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INTRODUCTION 


With the advent of the Space Transportation System and the era of career astro- 
nauts, who will spend appreciable portions of their careers in space, the problem of 
providing radiation protection increases in importance (ref. 1). For career astro- 
nauts and long-duration missions, the high-energy heavy-ion (HZE) component of galac- 
tic cosmic rays will become of major radiobiological significance, especially for 
nonregenerati ve tissues (ref. 2) • Since the range of these particles is large com- 
pared with a typical spacecraft wall thickness, nuclear attenuation and fragmentation 
appears to be a possible means of protection. 

In this work, initial estimates of fragmentation cross sections for a typical 

C 208^ C ra ^ nuc ^- eus ( 0 at 2.1 GeV/nucleon) colliding with a heavier target nucleus 

( Pb) are made and compared with experimental results (ref. 3). The fragmentation 
process is analyzed by using an abrasion— ablation collision model. The abrasion 
formalism utilized is described in detail in references 4 through 7. In order to 
complete our initial efforts at describing the fragmentation process, we adopted a 
simple ablation formalism where the excited projectile prefragment, which remains 
after the abrasion step, is treated as a compound nucleus and allowed to statis- 
tically decay by particle evaporation. The prefragment excitation energy is deter- 
mined from the geometric "clean cut" abrasion— ablation model of Bowman, Swiatecki, 
and Tsang (ref. 8). The compound nucleus particle evaporation probabilities are 
calculated by using the Monte Carlo program EVAP-4, developed at Oak Ridge National 
Laboratory (ref. 9). The resultant theoretical predictions are also compared with 
the predictions obtained from the geometric clean-cut model (ref. 8). 


PROJECTILE FRAGMENTATION 

The HZE particle fragmentation is assumed to take place in a two-stage formalism 
called the abrasion -ablation model. In this model, the projectile nuclei, moving at 
relativistic speeds, collide with stationary target nuclei. Those portions of their 
nuclear volumes which overlap are sheared away by the collision. This is the abra- 
sion step. The remaining projectile piece, called a prefragment, continues its tra- 
jectory with essentially its precollision velocity. As a result of the abrasion 
process, the prefragment is in an excited state and decays by the emission of gamma 
radiation and/or the evaporation of nuclear particles (ablation). The resultant 
isotope is the nuclear fragment whose cross section is measured. The abrasion part 
of the collision process is often analyzed from classical geometric considerations 
(refs. 8 and 10) or by using formal quantum scattering theory (refs. 4 to 7 and 11). 
The ablation part of the process may be analyzed by calculating the prefragment 
deexcitation through geometric arguments (ref. 8) or more sophisticated methods based 

upon Monte Carlo or intranuclear cascade techniques (refs. 10 through 13). 


Abrasion Cross Sections 

The cross sections a m for abrading m nucleons from the ^0 projectile were 
calculated by means of the quantum mechanical methods given in references 5 and 7. 
The appropriate oxygen and lead nuclear distribution parameters and nucleon -nucleon 
scattering parameters were taken from reference 14. The results are displayed in 



table I. (The symbols used in this paper are defined in a list after the refer- 
ences.) If n of the abraded nucleons are neutrons and z are protons, the 


TABLE I.- OPTICAL MODEL ABRASION CROSS SECTIONS FOR THE 
REACTION 16 0 + 208 Pb * m + x 

[Incident kinetic energy is 2.1 GeV/nucleon] 


Number of 
abraded 
nucleons, 
m 

Abrasion cross section, 

a , mb 
m 

1 

404 

2 

226 

3 

166 

4 

136 

5 

118 

6 

108 

7 

102 

8 

98 

9 

98 

10 

100 

1 1 

105 

12 

115 

1 3 

133 

14 

172 

15 

284 

16 

964 


resultant abrasion cross section, a nz can be obtained from a m by using a 
hyper geometric distribution 



where N is the total neutron number, Z is the total proton number, and the mass 
number is 


A = N + Z (2) 


It should be noted that the assumption of a hyper geometric distribution could result 
in inaccuracies, since it assumes there is no correlation at all between neutron and 
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proton distributions, and hence could lead to such unphysical results as abrading all 
neutrons or protons from a nucleus while leaving the remaining fragment intact. 


Prefragment Excitation Energies 

The excitation energy of the projectile prefragment following abrasion of m 
nucleons is calculated from the clean cut abrasion formalism of references 8 
and 10. For this model, the colliding nuclei are assumed to be uniform spheres of 
radii (i = P,T). In the collision, the overlapping volumes shear off so that 

the resultant projectile prefragment is a sphere with a cylindrical hole gouged out 
of it. The excitation energy is then determined by calculating the difference in 
surface area between the misshapen sphere and a perfect sphere of equal volume. This 
excess surface area A is given by (ref. 10) 


A = 4nRp[ 1 + P - (1 - F) 2 / 3 ] 


(3) 


where for 16 0 (R p = 3.32 fm) and 2US Pb (R^ = 7.04 fm), we have (ref. 10) 


208 


P = 


0.125(jiv) 1 /2^1 _ 


- 0.1 25 |o. 


W/ 2 (i_ 2 ) 


+ 1 


(M ! 


(4) 


and 


F = 0.75(1 - V )V»( 1 ^ - 0.125[3(1 - v) 1/2 - l] ^ 1 ~ PJ 


(5) 


wi th 


v 




( 6 ) 


P = 


b 


R + R 


(7) 


and 



(8) 
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Equations (4) and (5) are valid when the collision is peripheral (i.e., the two 
nuclear volumes do not completely overlap). For this case, the impact parameter b 
is restricted such that 


R - R < b < R + R 
TP TP 


( 9 ) 


If the collision is central, then the projectile nucleus volume completely over- 
laps the target nucleus volume (b < R T - R p ), and all the projectile nucleons are 
abraded. In this case, equations (4) and (5) are replaced by 


P = -1 


( 10 ) 


and 


F = 1 


(ID 


and there is no ablation of the projectile, since it was destroyed by the abrasion. 
If the excess surface area is used from equation (3), the excitation energy is 


E = A*E (12) 

exc s 

2 

where E g = 0.95 MeV/fm is the nuclear surface energy coefficient (refs. 8 and 10) 
from the liquid drop model of the nucleus. 


Ablation Factors 

Depending upon the excitation energy, the excited prefragment may decay by 
emitting one or more nucleons (protons or neutrons), composites (deuterons, tritons, 
^He, or alpha particles), or gamma rays. The probability a^j for forming a par- 
ticular fragment of type i as a result of the deexcitation of a prefragment of 
type j is obtained from the EVAP-4 computer code (ref. 9) by treating the prefrag- 
ment as a compound nucleus with an excitation energy given by equation (12). The 
final fragmentation cross section for production of the type i isotope is then 
given by 


a (z . ,A. ) = J2 a. . a . (z.,A.) 
F i i .11 abr i i 

j 


(13) 


The a __ (Z . , A . ) values are obtained from equation (1) by setting 
abr J J 


a (z . , A . ) = a 

abr D j nz 


( 14 ) 
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where 


z j 


= Z - z 


(15) 


and 


Aj = A - m 


(16) 


FRAGMENTATION RESULTS 

Table II displays the cross section results obtained from equations (1) and (13) 
for projectiles at 2.1 GeV/nucleon colliding with a 20 ®Pb target. Also displayed 
are the experimental isotope production cross sections from the fragmentation experi- 
ments described in reference 3. Despite the crude and unsophisticated nature of 


TABLE II.- ABRASION AND FRAGMENTATION CROSS SECTIONS FOR THE REACTION 

16 0 + 208 pfc| + A z + x 

[Incident kinetic energy is 2.1 GeV/nucleon] 


. A 

Species, Z 

a . (Z,A), mb 

abr 

a (Z , A) , mb 
F 

a (ref. 3), mb 

exp 

15 0 

202 

202 

135 ± 22 

14 0 

53 

0 

2.8 + 1.5 

13 0 

17 

17 


12 0 

5 

0 


n o 

1 .5 

0 


1 5 n 

202 

202 

202 + 26 

1 4 N 

120 

120 

71 + 22.5 

13 n 

66 

53 

17 ± 3 

12 n 

33 

0 


1 ’n 

15 

.6 


14 c 

53 

53 

12.3 ± 2.2 

13 c 

66 

0 

45.4 ± 8.3 

1 2q 

58 

132 

126+25 

"c 

43 

33 ! 

36.9 ± 5.7 

^o c 

26 

4.3 

7.21 ± 1 .40 

1 3 b 

17 

o 

0.7 + 0.4 

1 2 b 

33 

17 

3.98 ± 0.75 

1 L 

43 

33 1 

52.8 ± 5.9 

1 °B 

42 

85 

35.2 ± 11.3 

9 B 

35 

3 
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these fragmentation calculations/ the overall agreement between theory and experiment 
is quite good. Two main sources of error are the input excitation energy estimate 
and the assumption of a hypergeometric distribution. The excitation energy calcula- 
tion assumed that the excitation energy was independent of the charge or quantum 
states of the abraded nucleons and also ignored the surface diffuseness of the 
nuclei. The latter probably accounts for the overestimates of the 1 5 0 and 1 4 N cross 
sections, since the excitation energies of these species were too small for any par- 
ticle evaporation to occur. The most apparent example of the error due to the hyper- 
geometric distribution assumption is the overestimate of the 1 3 0 cross section. The 
use of this assumption resulted in a clearly unphysical 10-percent probability that 
the three abraded nucleons were all neutrons. This is also the likely source for the 
N overestimate, since 1 3 N is produced by evaporation of a proton from 14 0, which 
has a substantial cross section for removal of only two neutrons. The errors in the 
C and 1 C cross sections are attributable to the input excitation energy for 14 C. 
The C excitation energy (6 MeV) obtained from equation (12) is too small to ini- 
tiate particle evaporation. If it is increased, however, to «9 MeV, then neutron 
evaporation to 1 3 C can occur. This possibility is further supported by noting that 
the calculated a p for | 4 C is nearly equal to the sum of the experimental production 
cross sections for both 1 4 C and 13 C. 

Comparisons of the predictions of this work with those of the geometric 
abrasion-ablation model of reference 8 and the experimental results of reference 3 
are shewn in table III, where the cross sections for production of nitrogen, carbon. 


TABLE III.- ISOTOPE PRODUCTION CROSS SECTIONS FOR THE REACTION 


16 , 


>0 + 


208 


! Pb ■* Z + X 


[Incident kinetic energy is 2.1 GeV/nucleon] 


Process 

Isotope production cross sections, mb 

Geometric model (ref. 8) 

This work 

Experiment (ref. 3) 

0 ■* N 

419 

375 

290 + 35 

0 •> C 

286 

223 

228 ± 27 

0 -> B 

239 

138 

93 ± 13 


and boron isotopes from oxygen projectile nuclei are given. The production cross 
sections for a particular nuclear species (e.g., boron) were obtained by summing the 
contributions from all isotopes for that given species 


a p (z) =£c p (z,A) 


(17) 


As can be seen from table III, when the predictions of this work were compared with 
the esqperimental values of reference 3, closer agreement was obtained than by com- 
parison with predictions of the geometric model (ref. 8). 
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Future efforts at improving these and similar fragmentation calculations will 
center upon improving the method for estimating the prefragment excitation energy 
following the abrasion process. A more sophisticated quantum mechanical calculation 
using a weighted sum-rule method (ref. 1 1 ) is a possible improvement to this aspect 
of the theory. In addition, much of the neglected physics, including partitioning of 
the input excitation energy into rotational and internal degrees of freedom and cor- 
relations between neutron and proton distributions, needs to be incorporated. 


CONCLUDING REMARKS 

Predictions for the production of secondary nuclear species through the frag- 
menting of 16 0 projectile nuclei by 208 Pb target nuclei are presented and compared 
with experimental results. Although the theoretical calculations utilize very simple 
methods to treat the ablation effects, reasonable agreement with experimental data is 
noted. The major sources of disagreement appear to be attributable to the simple 
geometric model used to estimate the excitation energies of the projectile prefrag- 
ments and to the assumption of a hypergeometric distribution, which neglects correla- 
tion effects between the neutrons and protons in the nucleus. Clearly, much work 
remains to be done. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
December 12, 1983 
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SYMBOLS 


A 

b 

E 

exc 

E 

s 

F 

m 

N 


n 

P 

R P 

Rt 

z 


z 



a. . 

13 


P 


A 


a , (Z,A) 

abr 


exp 


m 

a 

nz 


nuclear mass number 

projectile impact parameter, fm 

projectile prefragment excitation energy, MeV 

nuclear surface energy coefficient, MeV 

defined in equations (5) and (11) 

number of abraded nucleons 

total number of nuclear neutrons 

number of abraded neutrons 

defined in equations (4) and (10) 

uniform nuclear radius of projectile, fm 

uniform nuclear radius of target, fm 

total number of nuclear protons 

number of abraded protons 

binomial coefficient 

probability of formation of type i fragment as a result of deexcitation 
of type j prefragment 

defined in equation (7) 

excess nuclear surface area, fm 2 

defined in equation (8) 

defined in equation (6) 

cross section for production of nucleus of type (Z,A) by abrasion, mb 
experimental cross section 
fragmentation cross section, mb 
cross section for abrading m nucleons, mb 
cross section for abrading n neutrons and z protons, mb 
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